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Abstract 

We successfully describe the HERA-data on diffractive deep inelastic scattering 
using a saturation model which has been applied in our earlier analysis of the in- 
clusive ep-scattering data. No further parameters are needed. Saturation already 
turned out to be essential in describing the transition from large to small values of 
Q 2 in inclusive scattering. It is even more important for diffractive processes and 
naturally leads to a constant ratio of the diffractive versus inclusive cross sections. 
We present an extensive discussion of our results as well as detailed comparison 
with data. 



1 Introduction 



In a recent analysis |1[ we introduced a model which provides a description of the transition 
between large and low Q 2 in inclusive lepton-proton deep inelastic scattering at low x. 
The idea behind our model is a phenomenon which we call a combined saturation at low 
Q 2 and low x. In this kinematical region the size of the virtual probe is of the order of 
the mean transverse distance between partons in the proton. The cross section for the 
interaction between the probe and the partons becomes large and multiple scattering has 
to be taken into account. These effects lead to saturation of the total cross section. We 
found that saturation occurs at low but still perturbative values of Q 2 (~ 1 — 2 GeV 2 
for x = 10~ 4 ). We therefore believe that saturation should be described by means of 
perturbative QCD (see also for example Ref. |2], [3|, [|, ||). 

The QCD-framework we use allows us to describe not only inclusive but diffractive 
processes as well. A general feature of diffraction is its strong sensitivity towards the 
infrared regime even for large Q 2 . The fact that diffraction has a strong soft component 
has already been noticed earlier, leading to the assumption that the Pomeron in diffraction 
ought to be soft. The idea of saturation, however, emphasizes the transition from hard to 
soft physics. As mentioned earlier saturation effects become already viable at rather hard 
scales and strongly suppress soft contributions in diffractive processes || . This mechanism 
leads to an effective enhancement of hard contributions and hence to an effective Pomeron 
intercept which lies above the original soft value. 

The important conclusion of this paper is that the concept of saturation leads to a 
good description of the diffractive data. Our approach has the important property that 
the inclusive and diffractive cross section have the same power-behavior in x. We have 
obtained these results without the use of any additional fitting parameters, i.e. we solely 
take the model as determined in Ref. [0 from the analysis of inclusive processes. The 
diffractive slope which we use is taken from the measurement at HERA. 

The plan of our paper is as follows. In Section 2 and 3 we recapitulate the results found 
in Ref. [||] and discuss qualitatively the basic features of saturation for both inclusive and 
diffractive processes. In Section 4 we provide the detailed formulae for our numerical 
analysis and discuss the comparison with the data from HI (7] and ZEUS || in Section 
5. In Section 6 we present the impact parameter version of our results and finish with 
concluding remarks in Section 7. The Appendix was added to enclose some details on 
the derivation of the cross section formulae for diffractive scattering. In particular it 
contains the computation of the relevant set of Feynman diagrams which contribute to 
the quark-antiquark-gluon final state. 

2 Saturation in inclusive processes 

Our starting point in Ref. |1J was the well established physical picture of small-x inter- 
actions in which the photon with virtuality Q 2 , emitted by a lepton, dissociates into a 
quark-antiquark pair far upstream of the nucleon (in the nucleon rest frame). The disso- 
ciation is then followed by the scattering of the quark-antiquark pair on the nucleon. In 
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this picture the relative transverse separation r of the qq pair and the longitudinal photon 
momentum fraction a of the quark (1 — a for the antiquark) are good degrees of freedom. 
In these variables the 7*p cross sections have the following factorized form ||, ID| 



ctt,l(x,Q 2 ) = Jd 2 v^da\^T,L(a,r)\ 2 a(x,r 2 ) , (1) 

where *&t,l is the squared photon wave function for the transverse (T) and longitudinally 
polarized (L) photons, given by 

|*T,L(a,r)| = T-^E e / (2) 

/ {4Q 2 a 2 (l-a) 2 K 2 (Qr). 

In the above formulae Kq and K\ are Mc Donald-Bessel functions and 

Q 2 = a (1 - a) Q 2 + m) . (3) 

The dynamics of saturation is embedded in the the effective dipole cross section a(x, r) 
which describes the interaction of the qq pair with a nucleon: 

a(x,r 2 ) =a |l-exp(- i ^- y )} , (4) 

where the x-dependent radius R is given by 

1 / x \ A / 2 

*o(*) = 77T7 -) • (5) 



GeV \XqJ 

The normalization o"o and the parameters x$ and A > of Ro(x) have been determined by 
a fit to all inclusive data on F 2 with x < 0.01 [p]]. (the detailed values of these parameters 
are quoted in Section 5). 

Saturation in the dipole cross section (0j) sets in when r ~ 2R Q , allowing a good 
description of the data at low Q 2 when 1/Q > R (see right plot in Fig. |l|). The detailed 
analysis of Eqs. (|T])-([|), presented in Ref. gives for small Q 2 

(TT ~ CTq . (6) 

For large Q 2 the dominant contribution comes from small dipole configurations with 
r ~ 2/Q <^ Ro (see left plot in Fig. |1|). In this case we have the usual situation of color 
transparency, a ~ r 2 , which gives the scaling behavior of the 7*p cross sections 

ot ~ 1/Q 2 • (7) 

More precise analysis leads to logarithmic in Q 2 modifications of the above estimations. To 
summarize, the saturation model naturally interpolates between the two different regimes 
of <7t described by Eqs. (j^) and (0). 
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Q 2 =10GeV 2 



Q 2 = 0.8 GeV 2 



■a 




Figure 1: The integrand of the inclusive cross section o~t in ([![) (solid lines) after the 
integration over a and the azimuthal angle, plotted for two values of Q 2 . The dotted lines 
show the dipole cross section ^). The dashed vertical lines correspond to the characteristic 
scales r = 2Rq and r = 2/Q. The values for 2R at a fixed energy W = 245 GeV are: 
0.36 fm for Q 2 = 10 GeV 2 and 0.25 fm for Q 2 = 0.8 GeV 2 . 

Saturation is characterized by a 'critical line'Q in the (x, Q 2 )-plane given by the equa- 
tion Q 2 = 1 / Rq(x) H], [| iJ. ft i s important to note that at very small x saturation effects 
become relevant at fairly high scales (Q 2 ~ 1 — 2 GeV 2 for HERA energies where one 
believes that perturbative QCD is valid. The critical line divides the phase space into two 
regions, the scaling region in which relation (0) is valid and the saturation region with 
the behavior given by Eq. ([]). 

The physical picture behind saturation is based on interpretation of the x-dependent 
radius Rq(x) as the mean separation of partons in the proton. We see from @ that when 
x decreases so does the mean seperation. Thus at low x the distribution of partons in the 
proton is no longer dilute when probed by a virtual photon of fixed resolution (~ 1/Q) 
and saturation sets in. This happens when the resolution of the probe equals to the mean 
separation, 1/Q = Rq(x), which condition defines the "critical line". As a result the dipol 
cross section becomes large and mulitple interactions become important. In other words, 
at low x proton appears to be black. The important result of our inclusive analysis is 
that blackening occurs already at rather short distances well below where 'soft dynamics' 
is supposed to set in, justifying the use of perturbative QCD. 

In the scaling region of large Q 2 the growth of the inclusive cross section is driven 
by the increase in the number of partons since the gluon density G(x) is proportional to 
1/-Rq(x) (see Section 4 for detailes of this relation). This grows is eventually tamed in our 
model by the mechanism of saturation. 



1 There is no phase transition or critical behavior present in our approach. 
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3 Saturation in diffractive deep inelastic scattering 



Inclusive j*p cross section at large Q 2 is dominated by the scaling region. Diffractive 
scattering on the other hand is essentially determined by the saturation region. In this 
case the dependence on x is controlled by the available phase space in the transverse 
momentum. This phase space grows proportional to 1/Rq(x) and leads to the same 
power behavior in x as was found for the inclusive cross section. It also means that the 
average transverse momentum of the diffractive final state will increase when x decreases. 
The process becomes 'harder' when x becomes smaller. Crucial for this picture to work 
is the scale invariance which in our approach is maintained by the lack of any additional 
cutoff on the transverse momenta of the final state. We will discuss this conclusion in 
more details below. 

In order to demonstrate the main features of saturation in diffraction we will confine 
our discussion in this section to the elastic scattering of the gg-pair as shown in Fig. 0a. 
Elastic gg-scattering dominates diffractive 7*p scattering for not to large values of the 
diffractive mass M. At large M, however, the emission of a gluon as depicted in Fig. |3|b 
becomes the dominant contribution. The cross section for elastic gg-scattering takes on 
the following form 



daS L 



dt 



t=o 



fd 2 r [da \y T , L (a,r)\ 2 & 

lOTT J JO 



2 / 2\ 

1 x, r 



with the same dipole cross section a as introduced for inclusive scattering. We account 
for the t— dependence by assuming an exponential dependence with the diffractive slope 
Bp- Thus the t-integrated diffractive cross section equals 



a D (x 



Q 2 ) = f° dte BDt 

J —oo 



da 



D 



dt 



t=0 



_L_ da D 



dt 



(9) 



t=o 



for both longitudinal and transverse photons. 

The distributions in r (qq dipole size) of the integrand for inclusive (Eq. ([I])) and 
diffractive (Eq. (§)) scattering at Q 2 = 10 GeV 2 are shown in Fig. |2|. The integrations over 
a and the azimuthal angle have been performed. The dotted lines denote the dipole cross 
section a. Comparing the two solid lines in Fig. 0a we see that for a typical inclusive event 
the main contribution is located around r ~ 2/Q <C 2Rq. The diffractive cross section 
on the other hand is dominated by the saturation region r ~ 2R . The importance of 
saturation for diffraction is illustrated in Fig. where we let a rise proportionally to 
r 2 . While this change has only little effect on the inclusive cross section, the diffractive 
cross section becomes strongly divergent One, in fact, needs an infrared cutoff - a new, 
additional scale R cut - to be introduce by hand. As a consequence an additional factor 
R 2 u1 -/Rq{x) emerges which leads a result reminiscent of the triple Regge approach where 
<j D ~ x~ 2X instead of a D ~ x~ x as we find with saturation. 

Fig. also illustrates the idea of Ref. that diffraction at small x is not a purely 
soft but semi-hard process. Let us assume for simplicity that the 'soft regime' begins at 
r = 4i?o- It becomes quite clear by comparing the two plots in Fig. |2| how strongly the 
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Q 2 =10GeV 2 



Q 2 =10GeV 2 




Figure 2: The integrands of the inclusive (Inc) and diffractive (DD) cross sections at 
Q 2 = 10 GeV 2 for the following two cases: (a) saturation according to Eq. (f^ (dotted 
line) and (b) no saturation, i.e. a ~ r 2 . 

soft contribution is suppressed due to saturation (blackness of the nucleon). The relative 
fraction of hard contributions (r < 2R ) is enhanced to almost 50%, making diffractive 
deep inelastic scattering a semi-hard process. A related issue is the smallness of the profile 
function in central collisions in pp scattering and its consequence for single diffraction. 

The following qualitative estimates will help to clarify the remarks about the impor- 
tance of saturation for diffraction. The wave function in Eq. (0) can be approximated by| 



(see also Ref. O, M) 



l*T,L(«,r)| 2 « ^Y.^{[^ + {l-^f]- 2 QHl-a)Q 2 r 2 <l] . (10) 
Att z j 1 r A 

The leading contribution is associated with the 'aligned jet' configuration. In the 7*- 
Pomeron CMS the scattering angle 9 is given by cos(#) = 1 — 2a, i.e. for a — > (1) we 
have 9 — > (it). The G-function in Eq. ([31]) enforces the condition that either a or 1 — a 
is smaller than l/(Q 2 r 2 ). We make use of this condition and the a «-> 1 — a symmetry to 
perform the a-integration in Eqs. (fj) and (||) and obtain 

o-{x,Q) « — — 2^e f — / —a(x,r) 

(11) 

a D (x,Q 2 ) « Eed /°° ^ a 2 (x, r 2 ) . 

v ' ^ ; 32n 2 B D y f Q 2 Ji/Q* r 4 v ' ; 

The lower limit is required, since the factor l/(Q 2 r 2 ) which results from the a-integration 
2 The relation Ki(x) ~ 1/x for x < 1 is used in Eq. (0) in the presented estimation. 
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should not exceed 1/4. We also approximate the dipole cross section (|j) by 

f ct r 2 /[AR 2 (x)} for r 2 < AR 2 (x) 
a » (12) 
[ cr for r 2 > 4i? 2 ,(a;) 

Inserting (|I2|) into ( [TT] ) gives after integration 



(13) 



16tt 2 5 d Y 7 4i^(i) Q 2 ' 



Thus we obtain an approximate constant ratio of the diffractve to inclusive cross sections 
similar to the exact result in Ref. [0 

°— « \ (14) 

a 8ttB d \n[R 2 (x) Q 2 } K ' 

If, on the other hand, we had used 

<^' g2) w ao 4^M (15) 

instead of (|TTD , i.e. no saturation, then a cutoff R 2 cut would be required leading to 

(16) 

~ 2.2**B D *f e * [ARl{x)] 2 Q 2 • 

The important point is that the inclusive cross section depends only weakly on i? cut 
whereas the diffractive cross section shows a strong dependence. We also realize that 
under the assumption (|T5| ) the diffractive cross section, being proportional to 1/Rq(x), 
rises at small x twice as strongly as the inclusive cross section (~ x~ 2X as mentioned 
earlier). The ratio (14) would be proportional to x~ x which is clearly not observed at 
HERA. 

To summarize, since the diffractive cross section is so sensitive to the infrared cutoff 
which is effectively given by 2Rq(x) one can conclude that diffraction directly probes the 
transition region. We will now turn to a full description of the diffractive deep inelastic 
scattering data from HERA. 
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4 Diffractive structure function in momentum space 
representation 



In this section we summarize the relevant contributions to the diffractive structure func- 
tion. We use the standard notation for the variables (3 = Q 2 /(M 2 + Q 2 ) and xjp = 
(M 2 + Q 2 )/(W 2 + Q 2 ) where M is the diffractive mass and W the total energy of the 
7*p-process. 

Before we start to compute the diffractive structure function it is useful to introduce 
the unintegrated gluon distribution T which is related to the effective dipole cross section 
(U) in the following way || |10[ : 



47T 

T 

4vr 2 



n 

11 



l-e"' 1 ] a s T(x,l 2 ) 
[l-JoM] a s F(x,l 2 t ) 



(17) 



A short calculation shows that with the following form for T 



3 o-Q 
4vr 2 



Rl(x 



p e -Ri{x)q 



(if 



one can indeed reproduce Eq. (||). At large Q 2 the usual gluon distribution G can be 
calculated from J 7 by a simple integration: 



xG(x,Q 2 ) = [ Q d£F(x,l 
Jo 



(TO 



3 (T 



47i 2 a s Rl(x)' 



1-(1 + Q 2 R 2 ) e~ Q2R o 
(Q 2 R%(x) » 1) 



(19) 



Important to note is the fact that at large Q 2 the gluon distribution exhibits a plain 
scaling behavior. The proper DGLAP evolution in Q 2 for the gluon can be added to 
our model, e.g. by treating relation (|i~9|) as the initial distribution for the linear DGLAP 
evolution equations. However the results of our model presented in Fig. 6 suggest that 
the Q 2 dependence of the data at low x values are properly accounted for in the presented 
approach and the additional gluonic evolution will only lead to a moderate improvement. 

We have three terms owing to the diffractive production of a quark-antiquark pair 
with transverse and longitudinally polarized photons and the emission of an extra gluon 
in the final state (Fig. ^j. The latter contribution is only known at present in certain 
approximations: strong ordering in the transverse momenta or strong ordering in the 
longitudinal momentum components. The first approximation is valid at a very large Q 2 
and finite diffractive masses, i.e. finite (3, and picks out the leading logarithm in Q 2 from 
the quark box. The second approximation is valid for large diffractive masses, i.e. small 
P, and finite Q 2 [0. Since the diffractive data range around (3 = 0.5 we will pursue the 
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Figure 3: Diffractive production of a qq-pair (left) and the emission of an additional gluon 
(right). 



first approximation and assume that the transverse momenta of the quarks compared to 
the gluon are much larger. 



For a detailed discussion of the derivation of the following formulae we refer to Ref. flcifl 
and only quote the result: 



xrFf^i&^xr) = ^- £ ej ^ £ da [a 2 + (1 



a 



(20) 



x 



dl 2 



1-2/5 + 



I 2 - (1 - 2(3) k 2 



[l 2 t + k 2 } 2 - 4(1 - p) 1} k 2 



and 



x P F° qq} (Q 2 7 f3,xip) 



-4r- t dak 2 (3 2 



(21) 



x 



r dlf 

y -pa s T(xip,l 



\J\ft + k 2 ] 2 - 4(1 - /3) Z t 2 P 

We have introduced the variable A; 2 which is defined as 

Q 2 _ fc 2 



a(l — a) 



(3 l-(3 



(22) 



and describes the mean virtuality of the exchange quark in the upper part of the diagram. 
Eq. (|22|) follows from the kinematics of the two-body final state. The variable a stems 



S 



from the Sudakov decomposition: k = aq' + (\k 2 \/2aq' ■ p)p + k t and q' = q + xp. The 
unintegrated structure function T is visualized in Fig. || as the lower blob. It is related to 
the inclusive F 2 by the optical theorem for zero momentum transfer t. In order to include 
the t-averaged distribution we have simply divided all expressions by the diffractive slope 
parameter Bp which has to be taken from the measurement, see Eq. (|9]). 

The third contribution takes on the formf]: 

9/3 ^ , 



%iP F $, qq - g }(Q 2 iPiXip) 



645 



x 



x 




<ft»^ln % 
2n \k 2 , 



dz 



(3 Z 2 (1 — Z) 



| J ~jT a sf(xiP, It) 



[(1 - 2z)k 2 - I 2 } 2 + 2z(l 



(23) 



z)k 4 



k 2 ^(l 2 + k 2 ) 2 - 4(1 



z) I 2 k 2 



In analogy to the previous formulae the variable k 2 expresses the mean virtuality of the 
exchanged gluon in the upper part of the right diagram (Fig. |]): 



k 



(24) 



The variable z represents the momentum fraction of the upper t-channel gluon with respect 
to the Pomeron momentum Xpp. It needs to be stressed that this formula was derived 
in the spirit of a leading log(Q 2 ) approximation which introduces uncertainties besides 
those related to the choice of a s . In this approximation the true kinematical constraints 
are not exactly fulfilled. The violation of these constraints, however, gives contributions 
which are sub-leading in the limit of very large Q 2 . An improvement can be achieved 
by an exact Monte Carlo integration. The exact treatment of the phase space, however, 
has to go along with the use of the exact matrix-element which is not known up to now. 
Similar analytic results can be found in Ref. ||14||. The main difference as compared to our 



approach is hidden in the unintegrated gluon distribution which in Ref. |I4[ is modeled 
by a heavy quark- ant iquark pair. 

There are two limits which are interesting to look at and which have been discussed 
in the literature: the first limit is the triple Regge limit (small f3) in which z can be set 
to zero in the square bracket of Eq. (EBI). This leads to 



xpF^ qqg} (Q 2 ,(3,xjp) 



x 



x 2 




2 ^_s_ 

2^ 



In 



'Q 2 



1 dz 

13 Z 2 



PY j [dl 2 t T( /2 , 



ml - k 2 ) + |e(A; 2 - / t 2 ) 



1\2 



'In Ref. |13| an overall factor of 2 was miscalculated and needed to be added. 



(25) 
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and agrees with results of Refs. [O, [15|, The other limit is I 2 <C k 2 and requires 



lower cutoff k 2 on k 2 : 



x P Fg qqg} (Q 2 ,P,x P ) 




d k 2 %L l n ( % 

kl 2tt ^ k 2 



dz 



p z 2 (1 — z) 



| a s xpG(x]p, k z 



(26) 



x 2 (1-zY (l + 2z) 



k 2 



This result and corresponding approximations for Eqs. (p0| ) and (pT|) have been derived 
earlier in Ref. ||17|| . They have been utilized in Ref. [TJj to perform a similar analysis of 
diffraction as presented in this paper. The approximations used in Ref. allow the direct 
implementation of the gluon structure function as given by standard parameterizations. 
The result is a too steep increase of the diffractive structure function with decreasing xjp. 
The exact formulae in conjunction with saturation give a much shallower behaviour which 
is in better agreement with the data (see below). 



5 Comparison with data 

Before we start our numerical investigation into diffractive scattering we would like to 
review the fit to the inclusive data The expression for the structure function Fi we 
have used in [|IJ was derived from Eq. (|I|) in combination with the saturation model (|J) 
quoted in Section 2. The parameters were found to be <To = 23.03 mb, A = 0.288 and 
Xq = 3.04- 10 -4 . These parameters enter into the diffractive cross section via the function 
T in Eqs. (p0|-|2"3|). To illustrate the quality of the fit we plot in Fig. | the structure 
function F 2 (x, Q 2 ) in different Q 2 bins in comparison with the data from HI |19| and 
ZEUS ]50] (see also for different comparison). 



The remaining integrations in Eqs. (pQ-|23|) have been performed numerically. We 
consider three light flavors and assume the diffractive slope parameter Bp = 6 GeV~ 2 
which is somewhat lower than the reported value of 7.1 GeV~ 2 plf . One has, however, to 



take into account some corrections due to double dissociation (dissociation of the proton) 
which can be roughly estimated by lowering the diffractive slope from 7.1 to 6 GeV~~ 2 . 
The coupling constant is kept fixed: a s = 0.2. 

Fig. [7| shows our result for the diffractive structure function xjpF D (xip, j3, Q 2 ) at fixed 
xp = 0.0042 plotted over (3 for various Q 2 together with data from ZEUS Fig. g 
contains similar plots with Hl-data for fixed xjp = 0.003 0. The three contributions 
( |2"Tp and (^3|) have been displayed separately in Fig. [7]. The important feature is the 



separation in three distinct regimes of small, medium and high (3 where the production 
of qqg, qq with transverse and qq with longitudinally polarized photons, respectively, is 
dominant. It was already argued in Ref. [|22| that this behavior is mainly due to the 



nature of the wave functions rather than the model we use. The relative strength of the 
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three contributions is fixed by QCD-color factors. The overall normalization, however, 
directly results from the saturation model without any fits to diffractive data. This fact 
is important to point out, since in Ref. [22| the overall and the relative normalization for 
the mentioned three contributions was fitted. One should note that there is no hard gluon 
component present in our approach (compare the analyses based on the concept of the 
'soft' Pomeron structure function |7], |28fl). 

The prediction of the xjp-dependence, besides the overall normalization, is an impor- 
tant consequence of the saturation model. In Fig. |9| and Fig. [XOj we compare our predictions 
with the data for xjpF D (x]p, j3, Q 2 ), now analyzed as a function of xjp for different values 
of (3 and Q 2 . Notice the good agreement, especially in the region of moderate and large 
values of (3 which corresponds to not too large values of the diffractive mass M. We 
also reproduce the change of the effective Pomeron intercept apuff) as a function of Q 2 
for different diffractive masses M, see Fig. |ll[ The effective intercept is related to the 
logarithmic xjp-slope n of X]pF D (xjp, /3,Q 2 ) through the relation: n = 1 — 2ajp( e //). At 
low masses M where the longitudinal part dominates the slope in xjp is slightly steeper 
due to the enhanced longitudinal part of the cross section. Using the effective Pomeron 
intercept means having incorporated shrinkage in the context of soft Regge phenomenol- 
ogy. The rise in Q 2 is again mainly caused by the longitudinal part. There is, however, 
another effect at work which lowers the intercept at small (3. The ggg-contribution has a 
logarithm ln(Q 2 /k 2 ) which is approximately equal to ln(Q 2 R (x) 2 ) . This term effectively 
lowers the intercept in the regime where qqg dominates, i.e. at small (3. 

In Fig. Owe show the ratio of the diffractive versus inclusive cross section as a function 



of W for different values of Q 2 and the diffractive mass M, in analogy to the analysis in 
Ref. . Thus for the presented analysis we have integrated Eqs. (|20D , ([H]) , (|23|) over the 
(3— values which correspond to the indicated ranges of M. The values of the inclusive cross 
section were taken from the analysis in Ref. . The ratio is almost constant over the entire 
range of Q 2 and W with a slight growth at small M caused by the longitudinal higher twist 
contribution. One can extract this behavior directly from the leading twist contributions 
of Eqs. (p0|) and (E3|) by simultaneous rescaling of the I 2 - and fc 2 -integration with respect 



to Rq. We have already discussed that the constant ratio is a particular feature of our 
saturation model and certainly deviates from the 'conventional' triple Regge approach. 



6 Diffractive structure function in impact parameter 
space representation 

We have started our discussion in impact parameter space because it provides a natural 
way to formulate saturation. For this reason we re-derive the formulae of Section 5 in 
impact parameter space. Moreover, the dipole formulation has its natural foundation in 
impact parameter space P3L A simple gg-pair represents an elementary color dipole 



which has an effective scattering cross section depending on the separation between the 
quark and antiquark. 



11 



We will briefly recall the wave function description for a gg-pair in impact parameter 
space using the conventions of Ref. [^] where the subscript (±, ±) denotes the photon- 
and quark helicity (complex notation): 

v/2 % C I I Y 

V'(+,+)(r,a) = a a{l ~ a)Q 2 a{\ - a)Q 2 r 2 ) - 

^(+,-)(r,a) = (1 - a) y/a(l - a)Q 2 a{\ - a)Q 2 r 2 ) (27) 

V(-,+)(r,a) = (1 - a) y/a(l - a)Q 2 a{\ - a)Q 2 r 2 ) j 



V2 



1 e 



■0(__)(r, a) = — ^— a ya(l - «)Q 2 -Ki(y a(l - a)Q 2 r 2 ) — . 
Ki is the MacDonald-Bessel function, and the variable r is conjugate to k t , i.e. 

V>(r,a) = j^e^^a) • ( 28 ) 
The longitudinal wave function reads: 



V(o,±)(r,a) = ^ a(l - a) Q K (ya(i - a)QV 2 ) . (29) 

The /9-integrated diffractive structure function can now be readily expressed in terms 
of the above wave function and the effective dipole cross section a [IlOL [TTJ : 



F /U g2 ' *) = E 4 / rf « [« 2 + a - «) 2 ] (3°) 



x a(l -a)Q 2 J d 2 r K 2 (^a(l - a)Q 2 r 2 ) a 2 (r } x) 



and 



3Q 2 f 1 

F{ lm} {Q\x) = 32n5BD E4 J Q daa(l-a) (31) 

x a(l -a)Q 2 J d 2 r K 2 (^Ja(l - a)Q 2 r 2 ) a 2 (r, x) . 

These two equation demonstrate the simplification one achieves in impact parameter 
space provided the distributions are totally integrated. They have already been quoted 
in Eq. (H) rewritten as diffractive cross section. The disadvantage, however, is that for 
differential distributions which depend on final state energies one has to transform back 
to momentum space as in the case of the /3-dependent structure function 

(3 2 f d 2 k t 

^ e (1-^11 02) 



1 _ 49 k ?- V 1 - PQ 

1 1-/3 Q 2 
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x J d 2 r fd 2 r' e ik ' (r - r,) a(r,x P ) a{r> \x P ) 

x 77 K > K > tfS**) ' 

We have made use of Eq. (p2|) to substitute a by 3 keeping k t fixed. For the longitudinally 
polarized photons we find 

xrF {l ^{Q ,/3,x F ) = £ e f ^— ^ J ^ k t 

k 2 /Q 2 ( A3 kf\ 

" 6, ( 1 -^t4 (33) 



4/3 Ay 2 V 1-/5Q 2 
1-/3 Q 2 

x /" d 2 r Jd 2 r' e ikt (r - r,) <r(r, a;,,) <7(r', x F ) 

This contribution is suppressed by an extra power in Q 2 and therefore is a higher twist con 
tribution. By using Eq. (|PTD one can directly transform Eqs. (|3~2]) and ([33|) into Eqs. (|20] 
and (1211). 



It should be noted that, when Eqs. ( |3"2"D and ( |3"3"D are integrated over 3, the argument 
xjp in a is simply substituted by x. This procedure is valid in the high energy approach as 
long as the dominant contribution is not concentrated at small 3. The /3-integration then 
leads from Eqs. ( j3"^) and ( |33"D back to Eqs. (|30| ) and ([!!]). In the case of a gluon in the 
final state one can no longer do a simple substitution but has to integrate the argument 
of a explicitly. 

We will discuss the impact parameterization of the ggg-final state in more detail. Our 
starting point is the wave function for the effective gluon dipole as described in [13[] (we 
use in this case the vector notation k 4 = {k\ , k 2 ) and m,n = 1,2 

/q(1 - a)Q- k 2 + a(l-a)Q 2 ^ 
1 k 2 t 5 mn - 2 k™k™ 



V^Q 2 kf + «q 2 

The second line of the previous equation is a consequence of the strong ordering condition 
which implies a < 1. The variable a has been introduced in analogy to Eq. (j22"D and is 
identical to zk 2 /(l — z)Q 2 , 

aQ * = = zk * , (35) 

where k t is the gluon transverse momentum in this case and k 2 describes the mean virtu- 
ally of the gluon in the upper t-channel. 
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The following relation illuminates the use of the wave function in momentum space. 
After the integration over the azimuth angle of l t one arrives at the core expression of 
Eq. (0) 

J ^ a s F(xjp, I 2 ) [2 ^ m >, k t ) - ^ m >, k t + l t ) - ^ mn (a, k t - L,)] (36) 



dl\ a. s T{x]p, I 



i 2 t V^Q 1 



dl 2 asT^xp^l 2 ) 



1 - 



+ 



2k 2 



I 2 aQ 2 



k 2 t + aQ 2 k 2 k 2 t 
[I 2 - k 2 + aQ 2 } 2 + 2k 2 aQ 2 



k 2 J[l 2 + k 2 + aQ 2 ] 2 -Al 2 k 2 



, 2 + (l-z) 2 + | 

[(1 - 2z)k 2 - I 2 ] 2 + 2z(l - z)k 4 



k 2 J(l 2 + k 2 ) 2 -A(l- z) l 2 t k 2 



\c m \c n ' 

5 mn - 2=^- 



The four terms ^ mn (a, k t )+ip mn (a, k t )—?p mn (a, k t + l t ) — ip mn (a, k t — l t ) represent the four 
possible ways of coupling the two t-channel gluons to the gluon dipole (without crossing 
in the i-channel). The Fourier transformation of the wave function leads to 



1 

2^ 



5 r 



aQ 2 K 2 (JaQ 2 r 2 ^ 



(37) 



Inserting the Fourier transform into the first line of Eq. ( |3lf) and using Eq. (|P7D we find 



J -^ta s F(x]p,l 2 t ) [2^ m >,k t ) - i/; mn (a,k t + k) - i> mn (a,k t -l t )} 



J d 2 r ^ mn («, r) e 4 ktr J ^ a s F(x F , I 2 ) (2 - e* '* r - e~ 4 ^ 



nil 



(38) 



2 / d 2 ri/; mn (a,r) e 



kt-r 



An 2 



a(xjp,r) . 



We can now rewrite Eq. (|23| ) in impact parameter space asQ 



a. r 1 dz 



1(3 z 



x 



d 2 k 



(2*) 



\i-z)Q' 

k 2 t 



e 



-9 + (f 

(1 - z)Q 2 - fc 2 



(1-*)S 

(39) 



J d 2 r J d 2 r' e ikt (r - r,) a(r, x F ) &(r', x F ) (V™ - 2^-) 



x <5 r 



r 



/2 



X 2 



1 - 2 



fc?r 2 ^2 



l-z 



2„r2 



kfr 



4 A missing factor 1/2 in the journal version was inserted in this update. 
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Again, a direct computation of Eq. fl39|) after substituting a according to Eq. (|T?D repro- 
duces the result of Eq. (^3|) . 

The impact parameter representation in Eqs. (|32| ) and (^) demonstrate the similarity 
of our approach and the semiclassical approach of Ref . [|(J . It suggests that the two-gluon 
exchange model can be extended to multi-gluon exchange without changing the basic 
analytic structure. The leading color tensors in the limit of large N c (number of colors) 
for a quark- and a gluon-loop with an arbitrary number of t-channel gluons attached to 
them are found to be identical up to an overall constant factor [27]. The large N c result 
differs only slightly from N c = 3 in the two-gluon exchange model and, hence, multi-gluon 
exchange is expected to give very similar results as the two-gluon exchange. 



7 Conclusions 

In our analysis we successfully describe diffractive deep inelastic scattering using the 
saturation model proposed in Ref. 0. This model reproduces quite accurately the f3- 
and xjp- distributions as measured by HI and ZEUS []7|, § without tuning or fitting any 
additional parameters. 

As demonstrated in Ref. |L| saturation naturally explains the transition of the inclu- 
sive structure function F 2 from high to low values of Q 2 . Diffractive scattering is even 
more effected by saturation (see Section 3). The constant ratio of the diffractive versus 
inclusive cross sections as observed at HERA is a direct consequence of saturation. It was 
also pointed out that soft contributions are significantly suppressed leading to a relative 
enhancement of semi-hard contributions. This fact allows the conclusion that diffraction 
in deep inelastic scattering is a semi-hard process ||. The effective Pomeron intercept is 
higher than expected from a 'soft' Pomeron approach |F| |28| . The /5-spectrum depends 
only weakly on the model and is therefore more universal. 

The model we choose for saturation is purely phenomenological. An alternative model 
without low-x saturation can be found in Ref. f31| . A completely theoretical framework 



involves non-linear QCD evolution equations as proposed in Refs. [0, [5], |5D[ . We believe, 
however, that our model represents the basic dynamics at very low x, since it allows us 
to describe a wide range of data in a satisfactory way. 

We can use our analysis to predict diffractive charm production. This requires the 
discussion of factorization, the introduction of diffractive parton distributions and the 
evolution of the diffractive final state. The detailed discussion of these topics will be 
presented elsewhere [f32|. 
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Appendix 



In this appendix we would like to recall the derivation of Eq. ( p3|) which represents the 
contribution due to the emission of an additional gluon |p9| . We choose light-cone gauge 
with the gauge fixing condition q' ■ A = (A is the gluon potential, q' = q + xp). The 
frame which naturally corresponds to this choice of gauge is the Breit frame, i.e. the 
frame in which the proton is fast moving. All quasi-Bremsstrahlungs gluons emitted from 
the gg-pair can be neglected. Those from the incoming partons on the other hand have 
to be taken into account. 

The polarization vector e for real gluons and the polarization tensor for the gluon 
propagator d Ufl read: 




Figure 4: Gluon radiation. 



Figure f| shows all the essential diagrams. The two diagrams to the left have a sim- 
ilar momentum structure and will be summed up right from the beginning whereas the 
diagram on the right will be calculated separately. The bottom line in all the diagrams 
represents a quark. It is accompanied by other 'spectator-quarks' which are not shown 
explicitly. The cut through the diagrams effectively subdivides the whole amplitude into 
two subprocesses. We will introduce effective three gluon couplings which are the sum of 
the original three gluon coupling and extra Bremsstrahlungs contributions (see Fig. [5]). 
These couplings and their analytic formulae represent the core of the whole calculation. 



17 



The blob at the top of the right t-channel gluon in Fig. |] indicates the simultaneous 
coupling of the t-channel gluon to the gg-pair which in color space combines into a gluon. 

Before starting the calculation one has to recall and make use of the kinematic assump- 
tions made in this approach. Firstly, there is the Regge limit with respect to the lower 
part of the diagram, i.e. the emitted gluon and the quark at the bottom have an invariant 
subenergy much larger than the diffractive mass M. The high energy assumption allows 
one to simplify the t-channel propagator as to 

dni + XrP) = r _ (l±^pi±^EEL (41) 

(A + %ip)p • q 

i p q '° 



(A + xf>)p ■ q' 



where the index p refers to the polarization at the upper end of the gluon line and o to 
the lower end. A corresponds to the Sudakov decomposition I = fyp + aiq' + l t where oti is 
fixed using the fact that the quark at the bottom is on-shell (ai ~ I 2 / s) . A itself is given 
through the on-shell condition of the intermediate s-channel gluon (l + xpp — k) 2 = and 
the final state gluon {xpp — k) 2 = 0: 



~ 1+ — 2// ' kf , . 

A = 1 — (42) 

a k s 



Oik 



(s = 2p-q'). Here the Sudakov representation of k enters with (3^ as free variable denoting 
the momentum fraction of the upper t-channel gluon with respect to the momentum p. 
Later on it will be substituted by z (z = [3k/ xp) which then denotes the momentum 
fraction of the t-channel gluon with respect to the Pomeron momentum. The contraction 
of q' a with the lower quark-gluon vertex gives roughly q' ■ p which cancels the same factor 
in the denominator of Eq. (f|1]). The remaining factor 1/(A + Xpp) in front of the vector 
l p t is large provided that xp is small. The other components of the polarization tensor 
d pa are negligible. All these properties are crucial in proving the /ct-factorization theorem. 
For the upper t-channel gluon the situation is different. In this case the corresponding 
tensor reads: 

d»°(k) = g>° - kPq \ + ^ (43) 

k ■ q' 

pa k p t q'° + q'Pkl + 2a k q' p q'a 

= 9t a ; 

PkV ■ q 

Due to the fact that the contraction of q' u downwards gives a factor x^p ■ q' which is not 
much larger than (3k, but of the same order, the term /3k in the denominator of Eq. (|43D 
is no longer enhanced as in Eq. (|4~T|) . However, a simplification is still possible, if one 
restricts oneself to the calculation of leading twist terms and keeps only the leading logs 
in Q 2 . Then, the transverse momenta of the quarks at the top of the diagram in Fig. ED 



and the gluon below are strongly ordered and all contributions with an extra inverse 
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power of the large quark transverse momentum are suppressed. This allows to set the 
transverse momentum k t along any of the quark lines to zero. Moreover, the projection 
of q' p with one of the upper quark-gluon vertices cancels or is sub-leading, and Eq. ([43]) 
may be reduced to: 



d pa (k) 



9t 



pa 



(44) 



This kind of technique is well known and has been applied in deriving the conventional 
Altarelli-Parisi splitting function. Therefore it is not surprising that the production of the 
gg-system is basically described by the AP-splitting function associated with the splitting 
of a gluon into two quarks accompanied by a logarithm in Q 2 jk\. Certainly, this approach 
is only valid for the transverse part of the cross section. The longitudinal part gives a 
next-to-leading log(Q 2 ) contribution which is not considered here. The coupling of the 
second gluon to the gg-system does not affect the dynamics within this system, but feels 
only the total color charge which is the same charge as carried by the first gluon. 











Figure 5: Effective triple gluon couplings. 



To summarize, the leading twist approach allows to factorize off the gg-system analo- 
gously to the conventional leading order DGLAP-scheme whereas in the lower part the kt- 
factorization theorem is applicable. All together, a local vertex may be extracted describ- 
ing the transition between the lower Pomeron exchange and the upper QCD-radiation. It 
is useful to rewrite Eq. (f44"l) in terms of transverse polarization vectors e t defined as 

9 r = (45) 

pol 

The sum has to be taken over the two helicity or polarization configurations in the trans- 
verse plane. d pa then reads: 

dF = £ P ( k ) £ a ( k ) (46) 

pol 
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with 

f( k ) = e° t - k J-^_ (47) 
PkP-q' K J 

Having in mind the previous discussion one can now start the calculation of the dia- 
grams in Fig. |j. The effective triple gluon vertex to the left of the first diagram gives the 
following contribution: 

2 p-e(k) l t ■ e(l + xjpp - k) 

- -5— e t ■ e(l + x P p - k) (48) 

- 2 p ■ e(l + x F p -k) If e t 

I 2 

- 2 p ■ e(l + xjpp - k) k t - e t 



The first three terms of Eq. ([4^) result from the ordinary three gluon coupling whereas 
the last is the sum of the two Bremsstrahlungs gluons as illustrated in the first row of 
Fig. [5| The momentum structure of these contributions is the same except the overall 
sign which is opposite. It is obvious that the two color tensors add up to the same tensor 
the ordinary three gluon coupling has. The overall color factor will be evaluated later. 
Here, only the correspondence between different color tensors is of interest but not the 
whole tensor itself. The right effective vertex in the first diagram of Fig. [| is different 
as it contains two s-channel gluons. Since these gluons are on-shell, the Ward identity 
l p A p = 0, where A p is the triple gluon coupling contracted with the gluon polarization 
vectors, may be used to change the t-channel polarization vector from It /Pi to p p . The 
resulting expression is: 

- 2 If e(xjpp -k) p- e(l + xjpp - k) 
+ 2 If e{l + xjpp — k) p-e(x]pp — k) 

— oikS e{l + xjpp — k) ■ e(xipp — k) (49) 



+ 2 p ■ e(Z + Xjpp — k) p-e(xjpp — k) 



I 2 

a k s 

Both pieces Eqs. (f|8D and (^9|) have to be combined and the sum over the transverse 
polarizations of the intermediate s-channel gluon has to be performed. The following 
equation will be used: 

£ 4(1 + xjpp - k) e»{l + x P p - k) = - gjT , (50) 

pol 

and products like p ■ e{k) will be reduced to —2k t ■ e t /Pk- Furthermore, the propagator 
1/k 2 = xjp/(a>ks) = (1 — z)/k 2 is introduced and Pi is expressed through Eq. (|42|) as well 
as the variable Pk is substituted by z (Pk = xjpz): 



2 kftt f Jf{k-h) 
xjp zk 2 \ (l t - k t ) 2 



I 2 

j^k t ■ e(xjpp - k) - If e(xjpp - k) 
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±U k - 

x P \ k 2 + % - 2l t ■ kt 



2 

+ — 



1 " z l t , 

If Ct H T2 fc * ' e * 

2 rC+ 



+ 2J* 



2 k t - e{x F p-k) 



(51) 



i (Zf - fcf) • £f 



Z 2 

■^fct • e(x F p - k) - If e(x]pp - k) 



k t ■ e{x P p - k) ' 
+ 2l t ■ e t p e t ■ t(x F p- k) 

I 2 

~Ti yokf e{x]pp- k) - If e{xjpp- k) 

{k - k t y Y k t 



+ 2 



k ■ (k ~ h) h ■ e(x]pp - k) (l t - k t ) ■ e(x]pp - k) 



(h - hY 



The next contribution has to be taken from the second diagram in Fig. [|. In this 
case the situation is slightly simpler compared to the first diagram, since only one effec- 
tive triple gluon vertex appears. Moreover, the upper t-channel gluon is attached to a 
quark line where the incoming and the outgoing quarks are on-shell with the consequence 
that the momentum of this gluon is purely transverse up to corrections proportional to 
the squared ratio of the gluon transverse momentum and the quark transverse momen- 
tum. This type of correction is sub-leading due to the strong ordering assumption. The 
polarization tensor simplifies in the following way: 



d pa (l + xpp - k) 



pfq" 7 
p ■ q' 



(52) 



The upper polarization vector was changed from — -J^^^a- to p p making use of the 
fact that the two quarks to the left and to the right are on-shell. In contrast to the 
first diagram in Figure |] the tensor gl" along the i-channel line gives only a sub-leading 
contributions due to the smallness of the longitudinal momentum. The special kinematic 
situation in the second diagram allows one to apply the eikonal approximation to the 
right quark-gluon vertex. The subsequent contraction with p p gives a factor which is 
cancelled by the residue of the 5-function corresponding to the intermediate quark, and 
the remaining factor is simply —1. The softness of the upper right t-channel gluon has no 
further dynamical effect except that the color charge of both quarks add up to the total 
color charge of the left t-channel gluon. Consequently, the color factor is identical to that 
of the first diagram in Figure [§. After all, one finds for this diagram: 



- ~r Iftt If ^(xjpp - k) (53) 
Pk 

2 I 2 

+ ~r h ■ e* V ■ A x tpV ~ k) 

Pk & k s 

Inserting the propagator l/(h — kt) 2 and substituting j3 k as well as a k one finally comes 
to: 

2 1 l t -e t _ 

k 2 t 



xtp z (l t - k t Y 



± k t ■ e(x]pp -k) - If e(x]pp - k) 



(54) 
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In the following step the two expressions (pl|) and (J5J) will be added and the result 
integrated over the azimuth angle between l t and k t . A lot of cancellations occur and the 
final expression is rather short: 

' ' [ , 2 + (l-,) 2 + | (55) 



2x P z{l-z) \ v ' k 2 

[(1 - 2z)k 2 - I 2 } 2 + 2^(1 - z)k 



k 2 J(k 2 + IX) 2 - 4(1 - z)l 2 t k 2 



e t ■ e t (x]pp - k) . 

Recalling the fact that only the amplitude has been considered, the calculation of the cross 
section requires to take the square of Expr. (|55|). In doing so one has to sum over the 
final state polarizations which leads to a contraction of the vector et with its conjugate. 
In the end the transverse part of the 7-matrices in the lower edges of the quark-box are 
contacted as well (see Fig. £|). 

Moving on to the final diagram (Fig. |]) one encounters a similar situation as in the 
case of the second diagram of the same figure. The right t-channel gluon is soft in the 
sense that its momentum is small compared to the quark momenta. It has no dynamical 
effect except that the color charge adds up as before, so that the final color factor is 
identical to that in the first two diagrams of Fig. [|. What remains is the calculation of 
the left effective triple gluon vertex. This has to be performed in a similar way as in the 
case of the left vertex in the second diagram: 

2 p ■ e(l + k) If e(x]pp — k) 

+ l l±2h^l e(l + k)-e(x F p-k) (56) 
Pi + xp 

— 2 p-e(xipp — k) l t - e(l + k) 

I 2 

- 2 p ■ t{xjpp - k) (l t + k t ) ■ e(l + k) 



(A + fa)ot k s 

The last term in Eq. (|56|) summarizes the contribution of the Bremsstrahlungs gluons 
associated with the effective triple gluon coupling. As was argued before the longitudinal 
momentum of the right soft t-channel gluon is negligible and fii equals zero. The momen- 
tum of the upper left t-channel gluon does not reduce to its transverse component, but 
includes the non-negligible longitudinal fraction z of the Pomeron momentum. Therefore, 
the propagator 1/(1 + k) 2 transforms into l/(zk 2 + (l t + k t ) 2 ). Introducing this propagator 
into Eq. (|56| ) and substituting fa = zxp as well as a^s = k 2 /xjp one finds: 

2 1 (l t + k t ) ■ e(l + k)lf t(xjpp - k) 



\2 



Xjp z zk 2 + (l t + kt) 

_2_ 1 lj_ (l t + k t ) ■ e(l + k)kf e(xjpp - k) 

Xtp z k 2 zk 2 + (l t + kt) 2 

2 1 h ■ e(l + k)kf e(x]pp - k) 

x F 1-z zk 2 + (l t + k t ) 2 [ } 

1 (h + k t ) 2 - k 2 

xp zk 2 + (l t + k t ) 2 
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Once more one has to integrate over the azimuth angle between l t and k t with the re- 
markable outcome that the resulting expression is identical to Eq. fl55|) : 



In other words, the sum of the first two diagrams in Fig. |4] is identical to the third diagram 
bearing in mind that the light cone gauge with the condition q' ■ A — was used. One 
should remind that the amplitude was calculated in the high energy asymptotic region 
where the real parts of the s-channel and w-channel contributions cancel due to the even 
signature of the color singlet exchange. (The w-channel contribution corresponds to the 
crossing of the two lower t-channel gluons in Fig. Hence, the imaginary part gives 
the leading part and was calculated taking the s-channel discontinuity, i.e. cutting the 
diagrams. However, the cut diagram gives twice the imaginary part and one has to divide 
the final result by 2. 

The structure in Eq. ( |58"D has been used in Eq. (^). The wave function in Eq. fl34|) 
cannot be extracted directly from the diagrams discussed here, but was constructed such 
that it reproduces the same results. 




(58) 



[(1 - 2z)k 2 - l 2 t ] 2 + 2z{l - z)k 4 



£t • e t (x]pp-k) 



k 2 ^(k 2 + l 2 ) 2 -4:(l-z)l 2 k 2 
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Figure 6: The results (solid lines) of the fit to the inclusive HERA data on F 2 for different 
values of Q 2 , using the model of [1] with saturation. 
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ZEUS 1994 




Figure 7: The diffractive structure function xpF D (x]p, ft, Q 2 ) for xp = 0.0042 as a func- 
tion of (3. The dashed lines show the qq contribution for transverse photons (|20|) , the 
dot-dashed lines correspond to the contribution from longitudinal photons fl2"T| ) and the 
dotted lines illustrate the qqg component (|23|). The solid line is the total contribution 
and the data are from ZEUS. 
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Figure 8: The same comparison as in Fig. |7]but with HI data. Only the total contribution 
is shown (solid lines). 
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Figure 11: The effective Pomeron slope as defined in the text as a function of Q 2 for two 
values of the diffractive mass M. 
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Figure 12: The ratio of the diffractive versus the inclusive cross sections as a function of 
W for different values of Q 2 and the diffractive mass Mx- 
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